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Abstract 



"^ We study a ferinionic Dark Matter particle carrying magnetic dipole moment and 

analyze its impact on direct detection experiments. In particular we show that it 
^ can accommodate the DAMA, CoGeNT and CRESST experimental results. As- 

suming conservative bounds, this candidate is shown not to be ruled out by the 
CDMS, XENON and PICASSO experiments. We offer an analytic understand- 
ing of how the long-range interaction modifies the experimental allowed regions, 
r-| in the cross section versus Dark Matter mass parameter space, with respect to 

Qh the typically assumed contact interaction. Finally, in the context of a symmetric 

Q^ Dark Matter sector, we determine the associated thermal relic density, and fur- 

ther provide relevant constraints imposed by indirect searches and colliders. 
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1 Introduction 

The presence of a relatively big dark component in the matter content of the Universe 
has been now established, and part of the scientific community is trying to unravel its 
still enigmatic properties. Despite the conspicuous number of experiments that probed the 
gravitational effects of this Dark Matter (DM) on the known matter, we still lack a more 
direct evidence of what it actually is and what are its properties. 

Under the assumption that the DM is constituted by one or more unknown particles, 
several experiments have attempted recently to probe its microscopic nature and its interac- 
tions. Different search strategies are possible. For indirect detection experiments, traces of 
DM annihilations to known particles are looked for in cosmic rays or in the spectrum of solar 
neutrinos. Collider searches are also performed, in case the DM can be produced in sizable 
amounts to have an impact in the channels featuring missing energy. Direct detection exper- 
iments rely instead on direct scatterings of DM particles off nuclei of ordinary matter, inside 
shielded detectors placed underground in order to diminish the cosmic rays background. 

DM direct detection experiments are providing exciting results. For example the 
DAMA/Nal collaboration [1] has claimed to have observed the expected annual modulation 
of the DM induced nuclear recoil rate [2,3], due to the rotation of the Earth around the Sun. 
The rotation, in fact, causes a different value of the flux of DM particles hitting the Earth 
depending on the different periods of the year. The upgraded DAMA/LIBRA detector has 
further confirmed [4] the earlier result adding much more statistics, and it has reached a 
significance of 8.9a C.L. for the cumulative exposure [5]. 

Interestingly the DAMA annual modulation effect has been shown to be compatible with 
a DM interpretation which, for the case of coherent spin-independent scattering, leads to a 
range of DM masses spanning from a few GeV up to a few hundred GeV, and cross sections 
between 10~^^ cm^ to 10~^^ cm^ [1,4,5], with some noticeable differences due to the galactic 
halo modeling [6,7]. 

More recently, the CoGeNT experiment first reported an irreducible excess in the count- 
ing rate [8], which could also be in principle ascribed to a DM signal. In the last months, the 
same experiment reported an additional analysis which shows that the time-series of their 
rate is actually compatible with an annual modulation effect [9]. The evidence of such a 
modulation for CoGeNT is at the level of 2.8a C.L. 

Also the CRESST collaboration observed an excess [10]. In particular, 67 counts were 
found in the DM acceptance region, where the estimated background from leakage of e/7 
events, neutrons, a particles and recoiling nuclei in a decays, is not sufficient to account for 
all the observed events. The analysis made by the collaboration rejects the background-only 
hypothesis at more than Aa [10]. 

The interesting feature is that the DAMA and CoGeNT results appear to be compatible 
for relatively light DM particles, in the few GeV to tens of GeV mass range and coherent 
scattering cross section around 10~^° cm^. CRESST points somehow to larger DM masses, 
but it is still compatible with the range determined by the other two experiments. The 
actual relevant range of masses and cross sections depends on assumptions of the galactic 
DM properties, namely the velocity distribution function and the local DM density [7] . Other 
relevant analyses can be found in Refs. [11-28]. 

The CDMS and XENONIOO experiments have recently reported a small number of events 
which pass all the selection cuts. Specifically, they have 2 events for CDMS [29] and 6 
events for XENONIOO, reduced to 3 events after post-selection analysis [30], which are 



still too few to be interpreted as potential DM signal. They can therefore provide upper 
bounds on the DM scattering cross section. These constraints seem to set severe bounds 
on the DM parameter space, and therefore to rule out much of the allowed regions of the 
other experiments. Very recently the PICASSO experiment [31] seems to provide even more 
stringent constraints for the low mass DM interpretation of DAMA, CoGeNT and CRESST. 

However, there are at least two caveats when interpreting the results from the experiments 
mentioned above. The first is that one has to pay attention to the fine details associated 
to the results quoted by each experiment, since a number of factors can affect the outcome. 
For example the actual response of the XENON and CDMS detectors can be uncertain and 
model dependent for a low energy signal [32-34]. Moreover, for crystal detectors like DAMA, 
an additional source of uncertainty is provided by the presence of an unknown fraction of 
nuclear recoils undergoing channeling. This effect is currently being investigated [35-37]. 
If confirmed, it would lead to a significant shift of the DAMA allowed regions in the DM 
parameter space. Another source of uncertainty is associated with the details of the nuclei 
form factors for each experiment. 

The second caveat is associated to the interpretation of the actual data within a 
very simple-minded model of the DM interaction with nuclei. For example, for the spin- 
independent case it is often assumed that DM couples via a contact interaction with equal 
strength to the proton and neutron. Upon relaxing some of these assumptions, less stringent 
constraints can be drawn. 

In this work we perform a thorough investigation of the effects of a magnetic dipole 
moment of the DM on the direct detection experiments. We refer to [23,38-50] for a limited 
sample of the earlier literature. 

A fermionic DM particle with a magnetic dipole moment can arise for example in models 
of composite DM [51-53], or more generally in models featuring a new strong interaction (see 
e.g. [54]). The application to cosmic rays as well as a complete discussion also of the light 
asymmetric DM paradigm, and last, but not the least, the first link to the study of grand 
unifications using DM appeared in [55]. Recently it was shown, for the first time, via first 
principle lattice simulations [56] that one has, indeed, for minimal models of Technicolor the 
correct pattern of chiral symmetry breaking leading to natural candidates of light composite 
DM. Despite the possibility to link the DM issue to other unresolved fundamental questions 
like e.g. the (dynamical) breaking of the Electroweak symmetry or to the grand unification 
problem, we investigate this DM candidate in a model independent way addressing only 
its DM phenomenology. The main feature of the present DM candidate is that it couples 
mainly to the photon and therefore interacts mostly with protons, and secondly and more 
importantly it has long-range type interactions modifying substantially the recoil spectrum 
with respect to a contact interaction. In fact, as pointed out in [57], low-energy threshold 
detectors like CoGeNT and DAMA are especially sensitive to this kind of interaction, and 
therefore the hope exists to get a better agreement among the various direct detection 
experiments. 

Here we will show how the DM magnetic moment interaction affects the interpretation 
of the different direct detection experiments. In Sec. 2 we present the DM interaction 
Lagrangian and provide some considerations concerning the scattering of DM off nuclei. 
This will allow us to determine the nuclear recoil rate for a given experiment. In Sec. 3 we 
show how the experimental favored regions and constraints are modified with respect to the 
naive contact interaction case. In Sec. 4 we describe the experimental data set we use and 



our statistical analysis performed to determine the favored regions and constraints for the 
DM candidate envisioned here. We summarize these results in Sec. 5. We determine the 
associated thermal relic density Sec. 6, and in Sec. 7 we report the constraints imposed by 
indirect searches, colliders and by the observations of compact stars. Finally, we conclude 
in Sec. 8. 



2 The event rate 

2.1 Kinematics 

When a DM particle scatters off a nucleus, depending on the DM properties, one can envision 
at least two distinct kinematics, the elastic and the inelastic. The elastic scattering is 
represented by 

X + U{A, Z)atrcst -> X + A^(^, ^)recoil , (1) 

while the inelastic is 

X + M{A, Z),trest ^X+ M{A, Z)recoil • (2) 

In (1) and (2), x and x! are two DM particle states, and A, Z are respectively the mass and 
atomic numbers of the nucleus M . In the detector rest frame, a DM particle with velocity v 
and mass m^ can scatter off a nucleus of mass rriNi causing it to recoil. The minimal velocity 
providing a recoil energy Er is: 



' rriN Er / /XxJV^ 

2/^2^ V T^N Er 



AEa)^d^^^[l + ^^], (3) 



where /x^at is the DM-nucleus reduced mass and 5 = m'^ — m^ is the mass splitting between 
X and x', and the equation above holds for 5 <^ m'^,m^. Elastic scattering occurs for 6 = 0, 
while 6^0 implies inelastic scattering. In this paper we will only consider the case of elastic 
scattering. 

2.2 Model and differential cross section 

We concentrate on the possibility of having a massless mediator of the DM-nuclei interaction 
able to yield long-range interactions absent in contact interactions usually assumed. The 
obvious candidate mediator in the SM is the photon while other exotic possibilities can be 
also envisioned such as a dark photon, see e.g. [57,58]. 

Assuming that the dark matter is either a neutral boson or fermion one can then use the 
language of the effective theories to select the relevant operators. Stability of the specific 
candidate can be easily ensured by charging it under an unbroken global symmetry. The 
most popular choices are either a new U{1) or a Z2 depending on the reality of the specific 
candidate DM field. 

In [59] one can find a general classification, up to dimension six, of all the possible 
interactions between a scalar DM and the SM fields. To this order one discovers that there 
is only one gauge-invariant operator coupling the DM complex boson to a photon, namely 



i{ip^X^)d,F'"' ; (4) 



DM-nucleus scattering mediated by this operator has been studied for instance in Ref. [18, 
60,61], where it is shown that it actually leads to a contact interaction and will henceforth 
not be used here. 

On the other hand, considering a fermionic DM x 0116 can show that the only gauge- 
invariant couplings to the photon, up to dimension five, are the magnetic and electric 
moment-mediated interactions: 

^M = -\k W.uXF^" , ^E =-\^d.x W.ul'xF^'' ■ (5) 

For a Majorana fermion both C^ and D^^ vanish identically, therefore x has to be a Dirac 
fermion. A^^ and d^ are the magnetic and electric dipole moment, respectively, and are usually 
expressed in units of e x cm. The scales Am = e/A^ and Ae = e/c?^ can be interpreted as 
the energy scales of the underlying interaction responsible for the associated operators to 
arise. For instance, in models in which the interactions (5) are explained by the DM being 
a bound state of charged particles, the A's could represent the compositeness scale. 
The differential cross sections for elastic scattering are [41] 
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where v is the speed of the DM particle in the Earth frame, a = e'^/Air ^ 1/137 is the 
fine structure constant, \p = e/2mp is the nuclear magneton and Fsi (-Fsd) denotes the 
spin-independent (spin-dependent) nuclear form factor which takes into account the finite 
dimension of the nucleus. Here 

\ isotopes / 

is the weighted dipole moment of the target [42], where /j, /Xj and Si are respectively the 
abundance, nuclear magnetic moment and spin of the i-th isotope; the values for these 
quantities are taken from [62,63] and agree exactly with the values provided in Fig. 1 of [42]. 
The differential cross section in (6) features both a spin- independent (SI) and a spin- 
dependent (SD) part. The SI part, neglecting for a moment the form factor, contains two 
terms: an energy dependent term with a E^^ drop-off of the cross section, and an energy in- 
dependent one. In contrast the common contact interactions only feature the constant term, 
typically with the 1/v^ dependence on the DM velocity; the SD part is of this kind. Notice 
that, for low enough energies (that might be also below threshold for some experiments, in 
principle), the interaction is always SI due to the Ej^^ divergence. As the recoil energy rises, 
the interaction becomes mostly SD for target nuclei with large magnetic moment, namely 
i9p^ 23]\]-g^ gj^^ i27j^ (Anuc/Ap = 4.55, 2.86, 3.33 respectively); the SD term is instead negligible 
for all the other nuclei. 



The differential cross section in (7), instead, has only an energy dependent term. We 
observe that the cross section in the electric case is enhanced with respect to the magnetic 
one by a factor 1/v'^ ~ 10^ translating in a value for d^ circa 10'^ times lower than the one 
for A-^, when one tries to fit the experiments. This is confirmed for instance by Ref. [41], 
where the authors find that, in order to fit the CoGeNT data alone, a Am of the order of the 
TeV is needed for a magnetic moment interaction, while for electric moment interaction Ae 
is around the PeV. Given that such a high scale is hardly reconcilable with other attempts 
to study DM and in general new physics, we will treat from now on only the magnetic dipole 
moment interaction O^. 

As for Fq,i and -Fsd, we use the nuclear form factors provided in Ref. [64]. For the SI 
interaction we have checked that Fq,i matches with the standard Helm form factor [65]. We 
recall that all the parameters used in the parameterization of the nuclear form factors may 
be affected by sizable uncertainties. 

2.3 Nuclear recoil rate 

The differential recoil rate of a detector can be defined as: 

dR ^^ f dauiv^En) 



where Nj- = N^/A is the total number of targets in the detector {Na is the Avogadro's 
number) and dn^ is the local number density of DM particles with velocities in the elemental 
volume d^v around v. This last factor can be expressed as a function of the DM velocity 
distribution f-E.{v) in the Earth frame, which is related to the DM velocity distribution in 
the galactic frame fciw) by the gahlean velocity transformation f-Eiy) = /g(^+^e(^)); here 
VE{t) is the time- dependent Earth (or detector) velocity with respect to the galactic frame. 
The prominent time-dependence (on the time-scale of an experiment) comes from the annual 
rotation of the Earth around the Sun, which is the origin of the annual modulation effect of 
the direct detection rate [2,3]. More specifically: 

Mt) = vg + vs + Mt) ■ (10) 

The galactic rotational velocity of our local system vq and the Sun's proper motion vs are 
basically aligned and their absolute values are vq = vq = 220 ± 50 km/s and vg = 12 km/s, 
while the Earth rotational velocity v^{t) has a size f® = 30 km/s, period of 1 year and 
phase such that it is aligned to vq around June 2^^*^ and it is inclined of an angle 7 ~ 60° 
with respect to the galactic plane. More details can be found, for instance, in Ref. [66]. 
Summarizing: 

dn^ = n^fE{v)d^v , (11) 

where n^ = ^^po/m^ is the local DM number density in the Galaxy and is determined by the 
local Dark Matter density po and, in general, by a scaling factor ^^^ which accounts for the 
possibility that the specific DM candidate under consideration does not represent the whole 
amount of DM. Here we assume ^^ = 1. In Eq. (11) the velocity distribution function needs 
to be properly normalized: this can be achieved by requiring that in the galactic frame 

d\fG{v) = l, (12) 



where fesc denotes the escape velocity of DM particles in the Milky Way. For definiteness, 
we will adopt here fesc = 650 km/s. 

When considering the differential cross section given in equation (6), the rate of nuclear 
recoils reduces to 
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with fmin(-E'R) given by Eq. (3). The detection rate is function of time through the velocity 
integrals X{vmm, t) and Ii(wmin, t) as a consequence of the annual motion of the Earth around 
the Sun. Their actual form depends on the velocity distribution function of the DM particles 
in the halo. In this paper we will consider an isothermal sphere density profile for the DM, 
whose velocity distribution function in the galactic frame is a truncated Maxwell-Boltzmann: 
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Under this assumption, and defining the normalized velocities tjq = {vg + 'ys)/'yo = Vq/vq, 
VE(t) = VE(t)/vo, r]ram{ER) = Vmm{ER)/vo and r]esc = Wcsc/^^0, the velocity integrals can be 
written analytically as [41] 
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where rj±{ER) = min(r7min(^i?) ±riE,Vesc)- 

Since in Eq. (10) the rotational velocity of the Earth around the Sun %, is relatively 
small compared to the main contribution represented by vq + vs, we can approximate VE(t) 
with its component directed toward the galactic center. We can then write [66] 



r]E{t) ~ ?70 + A?7 cos [2n{t - 0)/r] 



(20) 



where Ar/ = f0cos7/fo, with Ar/ ^ 7/0, and where cj) = 152.5 days (June 2°'^) is the phase 
and r = 365 days is the period of the Earth motion around the Sun. By means of Eq. (20) 
we can then expand the recoil rate, assuming that the velocity distribution is not strongly 
anisotropic: 
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m=ve 



To properly reproduce the recoil rate measured by the experiments, we should take into 
account the effect of partial recollection of the released energy (quenching), and the energy 
resolution of the detector: 



dR 



dE, 



det 



(Edet) = J dE' /C(Edet, E') J2^(^En = jy (22) 



Here the index i denotes different nuclear species in the detector, E^et is the detected energy 
and Qi are the quenching factors for each of the nuclear species. The function }C{Edet,E') 
reproduces the effect of the energy resolution of the detector; as is generally done, we assume 
for it a Gaussian behavior. 

Finally, the recoil rate of Eq. (22) must be averaged over the energy bins of the detector. 
For each energy bin k of width AEk we therefore define the unmodulated components of the 
rate Sok and the modulation amplitudes S^k as: 



Sok — . „ / dE^ct 



Smk — . „ / dE, 
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d dR 



AEk J AEk ^^^ d^' 
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At] . (24) 

m=v@ 

Sok and S^k are the relevant quantities that we use for the analysis of the experiments 
which address the annual modulation effect, namely DAMA and CoGeNT. For the other 
experiments, only the Sok are relevant. 

3 Theoretical predictions 

We attempt here an analytical study of the effect of the magnetic moment interaction on the 
observed differential rate. More precisely we provide a simple comparison between the results 
arising from the interaction studied here and the standard picture, i.e. the spin-independent 
coherent contact interaction, shown in Fig. 1. To this aim, Eq. (13) can be rewritten as 

^(t) = N^^^A^aXleiEn)XiVmUER),t)Fi,{En) , (25) 

where fi^n is the DM-nucleon reduced mass, aX^ plays the role of the spin-independent 
DM-proton cross section ap and we defined 

Q{En) = QsiiEa) f 1 + IttIt wIt) = ®si(^ii) (1 + r{En)) , (26) 
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Figure 1: DM-proton spin- independent interaction cross section o"p as a function of the Dark 
Matter mass m^, in the "standard" case of coherent contact interaction. The galactic halo 
has been assumed in the form of an isothermal sphere with velocity dispersion vq = 220 km/s 
and local density po = 0.3 GeV/cm'^. In this figure we show the allowed regions compatible 
with the annual modulation effects in DAMA and CoGeNT, as well as the region compatible 
with the CRESST excess, when interpreted as a DM signal. Specifically, the solid green 
contours denote the regions compatible with the DAMA annual modulation effect [4, 5] , 
in absence of channeling [35]. The short-dashed blue contour refers to the region derived 
from the CoGeNT annual modulation effect [9], when the bound from the unmodulated 
CoGeNT data is included. The dashed brown contours denote the regions compatible with 
the CRESST excess [10]. For all the data sets, the contours refer to regions where the 
absence of excess can be excluded with a C.L. of 7cr (outer region), 8a (inner region) for 
DAMA, la for CoGeNT and 3a (outer region), 4a (inner region) for CRESST. For XENON, 
the constraints refer to a threshold of 4 photoelectrons (published value [30] , lower line) and 
8 photoelectrons (our conservative estimate, upper line), as discussed in Sec. 4. The two 
lines for PICASSO enclose the uncertainty in the energy resolution [31]. 
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The function measures the deviation of the allowed regions and constraints with respect 
to the standard spin-independent picture, and the ratio r parametrizes the relevance of the 
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Figure 2: The function Q^En) (dashed hues), parametrizing the deviation of the allowed 
regions and constraints with respect to the standard case of contact interaction without 
isospin violation presented in Fig. 1. Each line is for one of the target elements used in 
the experiments taken into account in the text. The vertical lines indicate the low energy 
thresholds for the various experiments (apart from CoGeNT, for which we show the threshold 
of the second energy bin, where much of the signal is recorded). For DAMA two thresholds 
are shown, corresponding to the corrections due to the quenching factors for sodium (gNa = 
0.3, in the assumption of scattering mainly off Na) and iodine (gi = 0.09, for the assumption 
of scattering mainly off I). For XENON the 4 photoelectrons (4PE, published value [30]) 
and 8 photoelectrons (8PE, our conservative estimate) thresholds are shown, as discussed 
in Sec. 4. We have shown also the lowest threshold for the PICASSO experiment, although 
we refer to the main text for a more precise interpretation of its temperature dependent 
thresholds. The picture is reported for two values of the DM mass, m^ = 10 GeV (left 
panel) and m^ = 50 GeV (right panel). 

SD interaction respect to the SI one: if r > 1 the interaction is largely SD, while it is mostly 
SI for r < 1 (see Appendix A for a more detailed analysis). 

Since Qsi/^ depends itself on m-^, we expect an asymmetric shift of the favored regions 
and constraint lines in the (m^, cxp) plane due exclusively to the modified dynamics of the 
interaction considered here. There are only two free parameters in this model: m^, and Up (or 
alternatively A^). This is not true for instance in models with contact interaction in which 
the DM interacts differently with protons and with neutrons (the so-called "isospin violating" 
scenario, see e.g. [17, 18,22,58,67-74]). Indeed in such models the role of 6 is played by 
the expression {[Z + fn/ fp{A — Z)]/A) , and therefore there are three free parameters that 
must be fitted: cXp, m^, and the ratio of the DM-neutron and DM-proton couplings fn/fp- 

In Fig. 2 we show the behavior of G as a function of Eji both for small (left panel) and 
large (right panel) DM masses, considering several targets, Vq = 220 km/s and ve = Vq 
(being f0cos7 <^ ^0). In order to simplify the reading of the figures, we take the limit 
fcsc — ^ 00; one can figure out the effect of a finite escape velocity as a target-dependent 
cut-off at high recoil energies, that sets in when the minumum velocity required to have a 
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recoil above threshold exceeds fesc + ^e, the maximuni escape velocity in the Earth frame. 
There, the function G loses its meaning since the expected rate is zero both in the standard 
picture and in the case considered here. 

A first striking feature of is its overall magnitude: depending on the DM mass and 
on the experiment, we have in fact a suppression of the expected rate of roughly 1 to 4 
orders of magnitude. For nuclei in which the interaction is mostly SD (r > 1), provides a 
suppression of 1 to 3 orders of magnitude for the well known reason that the rate does not 
carry any A^ factor, contrarily to the usual SI case. Due to this fact, the SD part of the 
cross section is usually considered negligible respect to the SI one, but this turns out to be 
false in our case, for nuclei with large magnetic moment, since also the SI part is strongly 
suppressed. This suppression is due to the interplay of the last two factors in Eq. (27), 
while the first term plays only a marginal role given that A ~ 2Z roughly for all the target 
nuclei. The second term, 2^^^/mNEFt, enhances 9 by 2 to 4 orders of magnitude due to the 
presence of Eji, whose typical scale is few (tens of) keV and therefore very small compared 
to all the other mass scales involved. The last term Qq,i/X, containing the velocity integrals, 
gives instead a suppression of roughly 6 orders of magnitude. Given the measured rate, 
this overall suppression provided by the G function will reflect in the fit pointing to higher 
interaction cross sections, with respect to the standard case of Fig. 1. This fact will play a 
role in fitting the relic abundance (Sec. 6). 

The steep rise at low energies is always due to the long-range type of the SI part of 
the interaction, being the SD one of a contact type; in the particular case of magnetic 
dipole moment interaction considered here, the differential cross section diverges as 1/ En. 
At energies higher than the steep rise, for nuclei with large magnetic moment (r > 1), the 
interaction rapidly becomes of a contact type and therefore the function G exhibits a plateau 
whose value is ~ (AnucM^p)^. Instead, for the other nuclei for which r < 1, two regimes 
are possible, depending on the ratio fmin/'^E- These reflect the two trends assumed by the 
function^ C('^min) = '^livmin) /liy-aan) , namely a plateau C ~ 1.8fQ for t>min ^ v-p. and a rise 
C ~ "^min fo^ "^min ^ "i^E, with an intermediate value of C('^e) ~ 2.8 fg. 

Light DM particles require a higher minimum velocity to recoil compared to heavier ones. 
For the range of DM masses we are interested in, and for the range of recoil energies relevant 
for the direct DM search experiments, the value of f min/'^E is controlled by the ratio m^ jm-^. 
We will see that, depending on this ratio, the function Gsi can assume a constant behavior 
or scale as E'^ , while the function Gsd is always constant. 

3.1 Light Dark Matter 

For DM particles much lighter than the target nuclei, we have that v^^i > ve already at low 
energies for all the targets considered, and therefore 



gsi(^min(^fi),t) 
X(?7min(-Ei?),t) 



2 / r-. N r-. / 2miv + rn^ 



^min(-^ij) - Er 



2miym^ 



mNEn ^29) 



2m^ 



displaying the 



accordingly, G simplifies to (/ij(.„/(Amp)) Z"^ {nip/m^ + 2/3 (Anuc/Ap)' 

plateau shown in the left panel of Fig. 2. As expected, the plateau sets in earlier for heavy 

targets if the SI interaction dominates (Z^ {rrip/m^) ^ 2/3 (Anuc/Ap) ), while its starting 

^Here we ignore the time dependence, which gives only a negligible contribution. 
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point is independent on the target mass for SD interactions. We notice that, for high 
enough recoil energies, always gets to this constant behavior; it is interesting though that 
it exists a regime in which this happens at the low energies relevant for DM direct detection 
experiments. An indicative m^ ~ 10 GeV falls into this case, and as we will show in Sec. 5 
can fit the regions of the parameter space allowed by CoGeNT, DAMA-Na, CRESST-0 and 
CRESST-Cal 

For an experiment with a lower threshold above the steep rise, the differential rate does not 
feature any dependence on Er and it is therefore similar to the rate of spin-independent 
contact interactions, apart from a different functional dependence on m^. Notice that due 
to this different dependence on the mass of the DM, the rate is more suppressed for heavier 
DM particles, and therefore we expect the favored regions in the {m^, ap) plane as well as 
the exclusion lines to raise and slightly tilt at higher m^ (compared to the spin-independent 
contact interaction). 

Considering now that most of the signal in a given detector comes generally from the first 
energy bin^, i.e. close to the low energy threshold, we can roughly estimate, before doing 
any statistical analysis of the data, the shift of the allowed regions and constraints compared 
to the standard picture. For example, the expected rate of DM scatterings off Na nuclei is 
reduced by a factor ~ 10~^ with respect to the standard case (see Fig. 2), and therefore the 
DAMA favored region, when assumed that most of the signal comes from DM scattering 
on Na, is shifted up by ~ 10^ in the {m^, ap) plane. Taking now the DAMA-Na allowed 
region as benchmark, we can see how the other favorite regions and exclusion lines move with 
respect to it. CoGeNT moves a factor ~ 5 up, making the agreement with DAMA almost 
perfect. Concerning CRESST, the fit at small DM masses is due equally to O and Ca [10] 
in the standard scenario, from which we conclude that the cross section for DM magnetic 
dipole interaction with Ca is ~ 1.5 times bigger than that with O; the overall effect is an 
increase in ap, that lifts the favored region to better agree with DAMA and CoGeNT. In 
detail we expect that the CRESST allowed region moves a factor ~ 4 towards the DAMA-Na 
ballpark. Finally, CDMS and XENON move up by a factor ~ 6.5 and ~ 8.5 respectively, 
roughly irrespective of the choice of the threshold. This improves once again the agreement 
with the other experiments. For PICASSO the situation is different due to the special 
experimental setup in which the low energy threshold is a function of the temperature, and 
therefore one should be careful especially if the differential cross section is energy dependent. 
However, since the scattering in the PICASSO experiment is dominated by the dipole- 
dipole interaction which is of a contact type, a simple rescaling is again applicable, and we 
expect that the constraints will become a factor 2.5 more stringent (see left-panel of Fig. 2), 
compared to the standard picture. 

3.2 Heavy Dark Matter 

For higher DM masses, where for instance one can find the regions of the parameter space 
allowed by DAMA-I and CRESST-W, the effect of the long-range interaction turns to be 
very evident. As the minimal velocity for a given nuclear recoil v^i^ becomes smaller than 



^DAMA and CRESST are multi-target detectors and allowed regions at large DM mass correspond to 
scattering on I for DAMA and W for CRESST, while at small DM mass regions correspond to scattering on 
Na for DAMA and both O and Ca for CRESST. 

■^In CoGeNT, the larger part of the signal comes instead from the second energy bin. 
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We, the function C(^min) changes behavior. For fmin ^ ^e (i-e. m^ ^ uin] 
GsiiVminiE R),t) 



I{r]mm{ER),t) 



1.8.0^-^^''"^^ + "^- 



(30) 



In this regime, for the nuclei whose interaction is SI, scales as E]^^ for the whole range of 
the recoil energies relevant for the direct detection experiments. If En increases sufficiently 
enough, Wmin becomes eventually larger than ve and 6 assumes the constant behavior we 
described in the previous section. The difference in G between the low and high DM masses 
we depicted in Fig. (2) is that in the first case the plateau sets in at small recoil energies 
(smaller than the energy threshold of most experiments), whereas in the second case, the 
plateau sets in at large recoil energies (more than 25 keV). Instead for the other nuclei with 
large magnetic moment, the interaction is SD and as in the case of light DM the plateau 
manifests itself at low recoil energy. 

Due to the pronounced energy dependence in this case, the allowed regions and the 
constraints will shift in the {m^,ap) plane considerably more than in the case of low DM 
mass (still compared to the standard contact spin- independent interaction). By taking again 
the DAMA-Na allowed region as a benchmark from the right panel of Fig. (2), we estimate 
that both DAMA-I and CRESST-W shift more than one order of magnitude up in cXp. The 
same happens for CDMS and XENON. On the other hand PICASSO, like DAMA-Na, does 
not change its behavior respect to the light DM case: this is due to the fact that ^^F and 
^^Na enjoy mostly spin-dependent interactions, that depend only slightly on the DM mass. 

To summarize, taking into account the whole DM mass range, we see that, apart from an 
overall shift upwards in o"p, there is a general trend of the various experiments to "gather" , 
getting closer to each other, with respect to the standard case. This behavior is not homo- 
geneous in the DM mass and is more pronounced for masses above ~ 25 GeV. The effect 
is to favor a better fit compared to the standard case, with a better agreement of all the 
experiments. Even though this is true in general for any DM mass, we don't expect anyway 
to get a good agreement for higher masses: in that region in fact the XENON experiment 
rules out the signal featured by other experiments by several orders of magnitude in the 
standard picture (Fig. 1), and the shifts that take place in our case are not big enough to 
change this situation. 

Finally we comment on the dependence of the function G from the halo model considered, 
and in particular from the choice of the dispersion velocity vq. As we have seen, the transition 
between the two regimes discussed above is driven by the ratio Wmin/^E, and it is only 
important for the nuclei which experience SI interactions with a DM particle. Therefore 
in this case we expect that for larger values of ve ~ Vq, the plateau behavior of G sets 
in at higher recoil energies (or equivalently at smaller DM masses). This would make the 
shifts more pronounced even for light DM. In Sec. 5 and Figs. 3 and 4 we present the exact 
numerical results. 

4 Data sets and analysis technique 

In this section we discuss the techniques used to analyze the various data sets. In particular 
we adopt the approach of [57], and we summarize below the details of how we perform the 



^This is generally true provided we consider a recoil energy En within our interest: from a minimum 
threshold value of a few keV up to 25 — 30 keV. 
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fits to data and constraints from null results experiments. 

For DAMA, CoGeNT and CRESST, we test tlie null liypothesis (absence of signal on top 
of estimated background for CRESST and absence of modulation for DAMA and CoGeNT). 
From this we infer: 

i) the confidence level for the rejection of the null hypothesis (we find 8^ 9a for DAMA, 
1 ^ 2a for CoGeNT, and 4a for CRESST); 

ii) the domains in the relevant DM parameter space (defined by the DM mass m^ and the 
DM magnetic dipole moment A^) where the values of the likelihood function depart 
for more than na from the null hypothesis, and thus the corresponding evidence of the 
DM signal. We use n = 7, 8, n = 1, and n = 3, 4 for DAMA, CoGeNT, and CRESST, 
respectively [7]. 

Our statistical estimator is the likelihood function of detecting the observed number of 
events jC = ]^^ Ct, where the index i indicates the i-th energy bin in DAMA and CoGeNT, 
and the i-th detector in CRESST. For DAMA and CoGeNT £j are taken to follow a Gaussian 
distribution and for CRESST, since in this case the number of events in each sub-detector 
is low, a Poissonian one. Defining £bg as the likelihood of absence of signal, we assume the 
function y = —2 In C^g/C to be distributed as a x^-variable with one degree of freedom for a 
given value of the DM mass (notice that for DAMA and CoGeNT y reduces to y = xig ~ X^)- 
From the y function, the interval on A^ where the null hypothesis (i.e. A^ = 0) can be excluded 
at the chosen level of confidence are extracted: 7a (outer region) or 8a (inner region) for 
DAMA, la for CoGeNT and 3a (outer region), 4a (inner region) for CRESST. We then plot 
allowed regions in the {rn^, X^^) plane. 

We derive constraints from CDMS and XENONIOO with a similar likelihood function 
A = — 2 In C/ CTag] here C is the likelihood of detecting the observed number of events (2 for 
CDMS and 3 for XENONIOO), while in £bg the DM signal is not included. Both likelihoods 
are taken as Poissonian variables and A is assumed to follow a x^-distribution. For PICASSO 
we instead derive constraints by using a Ax^ method of the data points shown in fig. 5 of [31]. 
Bounds are conservatively shown at 5o" C.L. 

Concerning the experimental data sets and statistical methods, we refer to [57]; for 
completeness we summarize below for each experiment the most important ingredients used 
in the analysis. 

o DAMA: We use the entire set of DAMA/Nal [4] and DAMA/LIBRA [5] data, cor- 
responding to a cumulative exposure of 1.17 tonxyr. We analyze the modulation 
amplitudes S'^^ reported in Fig. 6 of Ref. [5] requiring that the DM contribution to 
the unmodulated component of the rate, 5*0, does not exceed the corresponding exper- 
imental value S'^'^ in the 2-^4 keV energy range. We compute 



( Q QCxp\2 / Q qcxpn2 

-21n£ = X^(A„ m,) = ^ ^^"^ ~ ^^^ > + ^^° ~ ^ ^ e{So - S^) , (31) 



^^ 



where ak and a are the experimental errors on S^^ and Sq"^^, respectively. The last 
term in Eq. (31) implements the upper bound on 5*0 by penalizing the likelihood when 
5*0 exceeds Sq^'^ with the Heaviside function 9. The detector energy resolution is 
parametrized by a Gaussian function of width aj.cs{E) = E{0A4:8/y/E + 0.0091) [75], 
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using for the quenching factor the central values quoted by the collaboration, namely 
Q'Na = 0.3 and gi = 0.09 [76]. We don't take into account the possibility of a nonzero 
channeling fraction 



o CoGeNT: We consider the time-series of the data, treating the measured total rate 
as a constraint. Similarly to the analysis done on the DAMA data, we define 



y = -21n£ = x^iK^m 



X) 



16 ( Q , , - C°^P ^ 16 /" Q ^ , _ Q'^'^P ^ 31 / c ccxpX 2 

2^ ~2 + Z^ —2 + 2^ —2 ^I'^Oj - ^Oi 

fe=l ^k k=l ^k j=i ^j 



(32) 



here S'mfc = l/At^ /^^ 5'mfc cos [27r(t — 0)/r] t/t, with /S.tk the temporal bin of the data, 
and 5^^ = i?^^ — (-R^fc)' where i?^^ is the total rate (taken from Fig. 4 of Ref. [9]) 
and (-R^fc ) is ^^^ annual average. The subscripts 1 and 2 in Eq. (32) refer to the first 
and second energy bins. The total rate in the 0.9 — 3.0 keVee energy-bin is computed 
by subtracting the rate in the 0.5 — 0.9 keVee bin to the rate in the 0.5 — 3.0 keVge bin, 
with a Gaussian propagation of the errors. S^'^ and cxj denote the experimental counts 
and the corresponding errors as given in Ref. [9] (31 energy bins in the interval 0.4 — 2 
keVee), after removal of the L-shell peaks but without removing any other background. 
The total fiducial mass is 330 g, the energy resolution is given by a Gaussian with width 
taken from [77], and the quenching factor below 10 keV is described by the relation 
^ = 0.2E]ji2 [78]. 

o CRESST: We compute the expected DM signal in each of the 8 CRESST detector 
modules. The acceptance regions and the number of observed events are provided 
in Table 1 of Ref. [10], and we derive background events according to estimates in 
Sec. 4 of Ref. [10]. A likelihood-ratio test yields a 4. la C.L. evidence for the best- 
fit of a DM signal over the background-only hypothesis, in good agreement with the 
result quoted by the collaboration. We use the published value of 730 kgxdays for the 
exposure and assume an even contribution among the different modules^ (each module 
accounts therefore for an exposure of 730/8 kgxdays); we consider moreover a constant 
efficiency. 

o CDMS: We use the "standard" 2009 CDMS-II resuhs based on Ge data [29]; these are 
obtained employing conservative nuclear recoil selection cuts and assuming an energy 
threshold of 10 keV. The total exposure is 612 kgxdays and we take the efficiency from 
the black curve of Fig. 5 in Ref. [79] with g ~ 1 as quenching factor.^ In spite of an 
expected background of 0.9 ± 0.2 events, two signal events were found in the 10 — 100 
keV energy interval [29] (we use these numbers to derive the constraints). 

o XENONIOO: We use the results presented in Ref. [30], with an exposure of 100.9 
days in a fiducial volume of 48 kg. After all the cuts, three events were reported in the 



^This is the same analysis performed in [57], although a value of 400/9 kg days for the exposure was 
erroneously reported in the text of the published version. 

^In the case of light DM, one can perform a similar analysis by using combined data from the CDMS and 
EDELWEISS experiments (see Fig. 1 in Ref. [80]), with basically the same results. 
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DM signal region in spite of an expected background of 1.8 ± 0.6 events. We model 
the data using a Poissonian distribution of photoelectrons, with a single-photoelectron 
resolution equal to 0.5. The shape of the Ces function is very crucial for such a low 
number of photoelectrons and for small DM masses. Following [57] we try to enclose 
a possible (but not exhaustive) uncertainty on the bounds derived for XENONIOO by 
adopting two different approaches: 

i) we adopt as threshold the published value of 4 photoelectrons and the nominal 
central value of C^s as shown in Fig. 1 of Ref . [30] , which relies heavily on linear 
extrapolation below 3 keVnr; 

ii) more conservatively, we raise the threshold for the photomultipliers to 8 photo- 
electrons: this value is the lowest one for which the analysis is nearly independent 
on the shape of C^s below 3 keVnr- 

Notice that these two approaches are not exhaustive of all the possible assumptions 
one can do to determine the XENONIOO response to light DM (for further discussion 
and considerations, see e.g. Refs. [32,33]). It appears therefore still preliminary, given 
these large uncertainties, to assume the bounds quoted by the collaboration as strictly 
firm. The 8 photoelectrons bound is less dependent on the Ccs extrapolation, and 
therefore, conservatively, we consider it as more appropriate. 

Finally, we follow Eqs. (13-16) in Ref. [81] to compute the expected signal. We derive 
upper bounds for both CDMS and XENON as mentioned at the beginning of this 
section. 

o PICASSO: The PICASSO experiment, located at SNOLAB [31], is very different from 
the ones discussed above; it is in fact based on the superheated droplet technique, a 
variant of the bubble chamber technique, to search for DM recoiling on ^^F in a C4F10 
target. The experimental procedure consists in measuring the acoustic signal released 
by the nucleation of a bubble as a function of the temperature T. Details of the 
detector principle can be found in [82,83]. Since bubble formation is only triggered 
above a certain energy threshold Eth(T), the spectrum of the particle- induced energy 
depositions can be constructed by varying the temperature. We compute the predicted 
DM rate as a function of Et^{T) from Eq. (3) of [31] and we compare such prediction 
with the experimental rate shown in Fig. 5 by using a Ax^ method. In this analysis 
we adopt two reference values of the parameter a(T) which describes the steepness of 
the energy threshold; namely we take a = (2.5, 7.5) in order to encapsulate as much 
as possible the experimental uncertainties [31]. We checked our result against the one 
of the collaboration given in Fig. 7 of Ref. [31], and we found an excellent agreement. 
Remaining coherent with our choice of being conservative we show here the result at 
5(7, as we did also for the other exclusion experiments. 

5 Fit to the direct detection experiments 

We analyze the direct detection data sets by using a standard isothermal halo model, which 
basically implies a truncated Maxwell-Boltzmann velocity distribution function (see discus- 
sion in Sec. 2.3). Since the response of direct detection experiments is quite sensitive to the 
DM distribution in the galactic halo [6] , especially in the scenario we are interested to study, 
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Figure 3: DM magnetic dipole moment A^ as a function of the Dark Matter mass m^. The 
galactic halo has been assumed in the form of an isothermal sphere with velocity dispersion 
Vq = 220 km/s and local density po = 0.3 GeV/cm^. Notations are the same as in Fig. 1; to 
match the two figures, one has to note that the role of ap is played here by «A^. The orange 
strip shows the values for (m^, A^) that fit the relic abundance Qum assuming a completely 
thermal DM production (see Sec. 6). 

we take into account uncertainties on the velocity dispersion vo, as discussed in Ref. [6]. We 
will use the three values Vq = 170, 220, 270 km/s, which bracket the uncertainty in the local 
rotational velocity. Let us notice that the value of the local DM density po is correlated to 
the adopted value of vq, as discussed e.g. in Ref. [6]. This corresponds to the model denoted 
as AO in [6] , and we adopt here the case of minimal halo, which implies lower values of the 
local DM density (since a fraction of the galactic potential is supported by the disk/bulge). 
In turn, this implies the adoption of po = 0.18,0.30,0.45 GeV/cm^ for vo = 170,220,270 
km/s, respectively [6]. 

Fig. 3 shows the constraints and the favored regions coming from DM direct detection 
experiments in the {m^, A^) plane. The galactic halo has been assumed to be in the form 
of an isothermal sphere with velocity dispersion vq = 220 km/s and local density po = 0.3 
GeV/cm^. The solid green contours denote the regions compatible with the DAMA annual 
modulation effect [4,5], in absence of channeling [35]. The short-dashed blue contour refers 
to the region derived from the CoGeNT annual modulation signal [9], when the bound from 
the unmodulated CoGeNT data is taken into account. The dashed brown contours denote 
the regions compatible with the CRESST excess [10]. For all the data sets, the contours refer 
to regions where the absence of modulation can be excluded with a C.L. of 7a (outer region), 
8(7 (inner region) for DAMA, Icr for CoGeNT, and the absence of an excess can be excluded 
at 3(T (outer region), Aa (inner region) for CRESST. Constraints derived by the null result 
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Figure 4: DM magnetic dipole moment A^ as a function of the Dark Matter mass m^. The 
galactic halo has been assumed in the form of an isothermal sphere with velocity dispersion 
vq = 170 km/s and local density po = 0.18 GeV/cm^ (left panel); vq = 270 km/s and local 
density po = 0.45 GeV/cm^ (right panel). Notations are the same as in Fig. 1; to match the 
two figures, one has to note that the role of Up is played here by aA^.The orange strip shows 
the values for {m^, X^) that fit the relic abundance I^dm) in the assumption of thermal DM 
production (see Sec. 6). 



experiments are shown at 5cr as gray, magenta and red dashed lines for CDMS, XENONIOO 
and PICASSO respectively. For the XENON detector, as discussed in the previous section, 
the constraints refer to thresholds of 4 and 8 photoelectrons, while for PICASSO we take 
into account the uncertainty in the intrinsic energy resolution of the employed detection 
technique [31]. 

We can see that, as expected from the discussion in Sec. 3, for a given value of m^, all 
the experiments determine a DM-proton cross section dp = aA^ that is higher than the one 
for the standard case of non isospin violating contact interaction depicted in Fig. 1; to give 
a benchmark value, DAMA points now to cXp ~ 1.5 x 10~^* cm^, about 10^ times above 
the standard case, as shown in Sec. 3. Moreover we observe the expected gathering of the 
various experiments, that get closer to each other with respect to the standard scenario, for 
an overall better agreement, and with DAMA, CoGeNT and CRESST featuring a larger 
overlap at low masses. Both the DAMA and CoGeNT regions point towards a DM mass in 
the 10 GeV ballpark (more specifically, from about 7 up to about 12 GeV) and DM magnetic 
dipole moment around 1.5 x 10~^^ ecm without exceeding the constraints, corresponding to 
an inverse mass energy scale of circa Am ~ 10 TeV. CRESST allows for even heavier DM 
masses, but still compatible with the range determined by the other two experiments. CDMS 
and XENON, on the other hand, exclude a smaller part of the parameter space with respect 
to the standard case. These two experiments do not exclude the overlapping region once one 
accepts our conservative choice of the XENON threshold at 8 photoelectrons. At 5a it is seen 
that PICASSO cannot exclude the common regions of the experiments reporting a signal. 
Finally, we also notice the occurrence of the expected flattening of the various experiments. 



meaning that the tilt featured in the standard case (higher mass regions pointing to lower 
(Tp) is very suppressed in the case of magnetic moment interaction, as explained in Sec. 3. 

These results depend on the galactic halo model assumed. The effect induced by the 
variation in the DM dispersion velocity is shown in the two panels of Fig. 4. In the case of 
fo = 170 km/s, the regions are not significantly modified as compared to the case Vq = 220 
km/s, except for the overall normalization due to the different values of the local DM density 
in the two cases. Larger dispersion velocities, instead, lead to the second regime discussed 
in Sec. 3 also at small DM masses (see right panel of Fig. 2). Due to the high energy 
dependence of the differential rate in this regime, experiments that use heavy targets and/or 
have high thresholds, get a suppression of the event rate, and point therefore to relatively 
higher values of the interaction cross section compared to the standard case. Since the DM 
signal is expected at lower recoil energies, smaller DM masses are now favorite. The right 
panel of Fig. 4 shows in fact that the allowed regions are shifted toward smaller DM masses 
for Vo = 270 km/s. 

6 Relic Abundance 

In addition to a satisfactory fitting of the direct detection data, the DM magnetic dipole 
interaction can accommodate the annihilation cross section appropriate for thermal produc- 
tion. Since it couples to the photon, the DM annihilates to any charged Standard Model 
particle- antiparticle pair, with mass lower than m^. Considering a DM particle with 10 
GeV mass, that as we argued fits well the direct detection experiments, the allowed primary 
channels are e~^e~ , fi'^fi^, t^t^ , uu, dd, ss, cc, bb and 77. 

The DM annihilation cross sections to Standard Model fermions and photons at tree level 
are 
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[\l{s\?>-l3'') + l2mls + ASm)ml)+dl{s\?>-l3^)-l2mls)] , (33) 
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where (3y_j = (1 — 4m^j/s)^/^, (3 = (3^(3f, Qf is the charge of the fermion and finally Nq = 1 
for leptons and A*"^ = 3 for quarks. In the non-relativistic limit, Eqs. (33) and (34) reduce 
to 
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where N = ^r Q'fNJj = 20/3 accounts for the number of degrees of freedom of the Standard 
Model fermions into which the DM can annihilate; the branching ratios BR(/) are defined 
as BR(/) = Q'JNIj/N. As one can see, for the low values of A^^m^ pointed by the direct 
detection experiments, the two photons final state is suppressed with respect to the fermionic 
one. Using a DM energy density ^dm^^ = 0.1126 ± 0.0036 measured by WMAP [84], we 
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show^ with an orange strip in the {m^, X^) plane of Figs. 3 and 4 the phase space where a 
thermal production of the magnetic DM is possible. 

It can be seen in Figs. 3 and 4 that the agreement between the allowed regions and 
the relic abundance depends strongly on the dispersion velocity Vq. The best agreement 
takes place for lower values of Vq. This is due to the fact that, as commented already in 
Sec. 5, po is linked to the dispersion velocity vq- By varying vq, po is modified accordingly, 
changing therefore the expected rate and the fitting values for the interaction cross sections, 
or in other words A;^. As a consequence, the favored regions and constraint lines of direct 
detection experiments move along the vertical axis in the {m^, A^) plane. For lower velocity 
dispersion, both the local DM density and the expected rate decrease. Therefore, the favored 
regions will point to a higher DM magnetic moment, getting closer to the relic abundance 
strip. The value of the relic abundance is obviously independent on the local DM density. 

The fact that the magnetic DM (with only two parameters) provides a good fit to the 
direct detection experiments and simultaneously can accommodate a thermal annihilation 
cross section (which is several orders of magnitude stronger than the DM-nucleon cross 
section in the context of "standard" contact interaction) is not easily met by other DM can- 
didates. Although there are candidates that can have a DM-nucleon cross section much lower 
than their thermal annihilation cross section, these candidates have usually spin-dependent 
interactions and cannot fit nicely the direct search experiments with a positive DM signal. 
For a typical candidate with spin-independent contact interactions, and for a strength of 
interaction that leads to a DM-nucleon cross section pointed by the direct detection exper- 
iments featuring a signal, the annihilation cross section is way too small to produce this 
candidate thermally. In order to match, a suppression mechanism for the direct detection 
event rate should take place. In the case of magnetic moment DM, this mechanism exists 
and has two different reasons for the SI and the SD parts of the interaction. As already 
explained in Sec. 3, this suppression is encoded into the function 9 in Eq. (26), that spans 
roughly from 10~^ to 10~^ for a 10 GeV DM (see Fig. 2). This suppression is a result between 
two competing factors for the SI case: the enhancement provided by the E]^^ dependence 
of the differential cross section, and the suppression provided by the kinetic integral; in the 
SD case, the suppression is instead due to the lack of the A^ enhancement usually present 
in the standard SI case. 

Notice that in the case of electric dipole DM one faces a very different situation: here the 
differential cross section has the same dependence on the velocity as in a contact interaction 
(Eq. (7)), and therefore there is no kinetic suppression, while on the contrary the rate 
is enhanced by the E]^^ dependence. Moreover, the annihilation to fermions is a p-wave 
process, Eq. (33), and therefore the annihilation cross section is suppressed, making the 
value of d^ needed to fit the relic abundance even bigger. 

Furthermore, it is worth to point out that neither asymmetric/mixed DM [86], nor os- 
cillating DM [87-89] can improve the agreement between the relic density and the allowed 
regions of direct DM searches in the context of magnetic DM interaction. 



'The figures are produced by solving numerically the Boltzmann equation 

Titot + SiJntot = -X (f^annl'rcl) ("-Lt " "cq) : (37) 

where ntot is the total number density of particles and antiparticles, H is the Hubble parameter and for 
(canni'rei) we usc the morc precise formula given in Eq. (3.8) of [85]. 
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Figure 5: DM magnetic dipole moment A^ as a function of the Dark Matter mass m^. The 
galactic halo has been assumed in the form of an isothermal sphere with velocity dispersion 
Vq = 220 km/s and local density po = 0.3 GeV/cm^. Notations are the same as in Fig. 1. The 
orange strip shows the values for {m^, A^) that fit the relic abundance I^dm, in the assumption 
of thermal DM production (see Sec. 6). The shaded regions refer to the constraints from 
7-ray lines, CMB, galactic 7-rays and dwarf galaxies. These constraints are only valid in the 
assumption of symmetric DM. The galactic photons constraint enforces total annihilation of 
the DM into bb, and therefore a less stringent constraint is expected for magnetic moment 
DM. 

7 Constraints 

Given the qualitative agreement between the direct detection allowed regions and the fit 
of the relic abundance, it is natural to ask whether constraints coming from other searches 
could limit the parameter space of the magnetic moment DM. This is even more pressing 
since the bounds coming from indirect DM searches are on the verge of constraining the 
thermal annihilation rate for light DM particles. We will discuss the constraints imposed by 
indirect searches, colliders and by the observations of compact stars, and then identify the 
most stringent ones. 

7.1 Epoch of reionization and CMB 

Strong constraints are imposed on DM annihilations from considering the effect on the 
generation of the CMB anisotropics at the epoch of recombination (at redshift ~ 1100) and 
their subsequent evolution down to the epoch of reionization. The actual physical effect 
of energy injection around the recombination epoch results in an increased amount of free 
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electrons, which survive to lower redshifts and affect the CMB anisotropics [90-93] . Detailed 
constraints have been recently derived in [94,95], based on the WMAP (7- year) and Atacama 
Cosmology Telescope 2008 data. The constraints are somewhat sensitive to the dominant 
DM annihilation channel: annihilation modes for which a portion of the energy is carried 
away by neutrinos or stored in protons have a smaller impact on the CMB; on the contrary 
the annihilation mode which produces directly e~^e~ is the most effective one. Usually the 
approach here is to consider 100% annihilation rate in a single final state; anyway, in our case 
several annihilation channels are open, and therefore we expect a smaller energy injection 
in the interstellar medium with respect to the case of annihilation only into e^e~ . In Fig. 5 
we reproduce the constraints as obtained in Refs. [93-95] considering now all the channels, 
each with its branching ratio defined in Sec. 6. 

7.2 Present epoch 7-rays 

For most of the DM annihilation modes, another relevant constraint is in fact imposed 
by the indirect DM searches in the present epoch. The DM constraints provided by the 
FERMI-LAT 7-ray data are particularly relevant as they are now cutting into the thermal 
annihilation cross section for low DM masses (< 30 GeV) and a variety of channels. 

In particular, dwarf satellite galaxies of the Milky Way are among the most promising 
targets for Dark Matter searches in 7-rays because of their large dynamical mass to light 
ratio and small expected background from astrophysical sources. No dwarf galaxy has been 
detected in 7-rays so far and stringent upper limits are placed on DM annihilation by applying 
a joint likelihood analysis to 10 satellite galaxies with 2 years of FERMI-LAT data, taking 
into account the uncertainty in the Dark Matter distribution in the satellites [96]. The limits 
are particularly strong for hadronic annihilation channels, and somewhat weaker for leptonic 
channels as diffusion of leptons out of these systems is poorly constrained. In our case, 
having both hadronic as well as leptonic annihilation channels, we expect again a weaker 
constraint with respect to the pure hadronic annihilation (see e.g. figure 2 in Ref. [96]). 

Other strong limits on annihilation channels are set by, for example, the 7-ray diffuse 
emission measurement at intermediate latitudes, which probes DM annihilation in our Milky 
Way halo [97-99]. In particular, the most recent limits come from 2 years of the FERMI- 
LAT data in the 5° ^ 6 ^ 15°, -80° ^ £ ^ 80° region [99], where b and i are the galactic 
latitude and longitude. Since bounds on all annihilation channels are not available with the 
latest data, we report only the most stringent one, coming from bb. 

The last kind of constraints that we can set on DM annihilation in the case of magnetic 
moment interaction comes from 7-ray lines. Indeed, given the possibility of annihilation in 
two photons, we expect a line in the cosmic photon spectrum at energies equal to the mass 
of the DM. Constraints on the annihilation cross section into photons can be drawn from 
the latest FERMI-LAT data [100]. For instance, the FERMI-LAT collaboration excludes at 
95% the annihilation {a^^^fj frci) < 0.5 x 10"^'^ cm^/s for a DM mass m^ = 30 GeV in 
the case of an isothermal halo. In our case we use the latest (preliminary) data on the flux 
from spectral lines (shown in slide 39 of [101]) to constrain the annihilation cross section of 
Eq. (36)^. We consider an all-sky region with the Galactic plane removed {\b\ > 10°), plus 
a 20° X 20° square region centered on the Galactic center [100]. Notice that, for small DM 
masses and dipole moments, the one-loop contribution becomes comparable with the tree 

*A similar study has been published in Ref. [102] for DM masses above 30 GeV. 
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level one; a rough estimate of the loop cross section in the non-relativistic limit is a^A^/47r, 
and therefore we expect it to become sizeable for Xyjn^ < (T'I'^ ~ 6 x 10~^. In the part of 
the parameter space probed by 7-ray lines searches the loop correction is negligible. 

We superimpose all these constraints in Fig. 5. We see that, apart from the 7 lines 
bound, they are somewhat stronger than the CMB one considered above. We keep however 
the latter as it is less model dependent. 

As shown in Fig. 4, a change in the DM local density modifies the direct searches results. 
This does not happen for indirect searches, as they are not very sensitive to po- Therefore 
the constraints shown in Fig. 5 also apply, unchanged, to this case. For lower DM dispersion 
velocity fo these constraints start to play an important role in cutting the parameter space 
for direct detection searches. 

7.3 Collider and other Astrophysical constraints 

The constraints we have addressed above apply to the case of symmetric (thermally pro- 
duced) DM. These constraints are not relevant in the case where DM is of asymmetric 
nature. There are two extra type of constraints: constraints imposed by collider searches, 
and constraints imposed by observations of compact stars such as white dwarfs and neutron 
stars. The collider constraints are applicable wether DM is symmetric or not, whereas the 
compact star constraints are valid only for asymmetric DM. 

The collider constraints emerge from the fact that for a given A^ and m^ that fit the 
direct DM searches experimental data (and even provide the proper DM annihilation cross 
section for thermal production), the cross section for production of a pair of DM particles 
in colliders is completely fixed. The processes that lead to constraints are mono-photon 
production (from initial or final state) plus missing energy due to the pair of DM particles 
in e^e~ collisions at LEP, or mono-jet production plus missing energy in proton-antiproton 
collision in Tevatron, or proton-proton collisions at LHC In the case of magnetic DM these 
constraints have been studied in [48,49] where it is found that the upper bound on A^ is 
safely above the range of values of A-^ relevant for the direct search experiments (for the 
range of m^ ~ 10 GeV). 

In the case of asymmetric DM, constraints can be imposed by compact star observa- 
tions based on the fact that a substantial number of captured DM particles might lead to 
gravitational collapse and formation of a black hole that can destroy the host star. The 
magnetic DM being a fermion evades the severe constraints on asymmetric bosonic DM 
based on neutron stars [103-105]. The constraints on self-DM cross section with Yukawa 
interactions presented in [106] are avoided for several reasons: firstly DM-DM interactions 
scale as fn'^X^ (leading to a typical DM-DM cross section of ~ 10~^^ cm^ which is much 
smaller than the constraint). Secondly the constraints are not directly applicable because 
the mediator of the magnetic DM is a massless photon (and not a massive mediator nec- 
essary for the constraint) and the DM-DM interaction is repulsive. The latter adds up to 
the effect of the Fermi pressure of the DM particles and therefore the amount of particles 
needed for gravitational collapse cannot be accumulated within the age of the Universe. A 
potential attractive photon interaction takes place in the symmetric case, which will however 
lead also to annihilation of the DM population inside the neutron star invalidating thus the 
constraint derived from black hole formation. 
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Finally the magnetic DM evades constraints on the spin-dependent part of the cross sec- 
tion imposed by observations of white dwarfs [107]. Although these constraints are typically 
weaker than the ones derived from direct searches at the range m^ ~ 10 GeV, they could 
become stricter if spin-dependent interactions scale as some (positive) power of the recoil 
energy. Since DM particles acquire high velocities when entering the white dwarf, such con- 
straints could in principle exclude such a candidate. However, as it was pointed out in [41], 
the spin-dependent cross section does not scale with the recoil energy and therefore the white 
dwarf constraints can be safely ignored. 

8 Conclusions 

We investigated a fermionic Dark Matter particle carrying magnetic dipole moment and 
analyzed its impact on direct detection experiments. We provided an analytic understanding 
of how the photon induced long-range interaction neatly modifies the experimental allowed 
regions with respect to the typically assumed contact interaction. We showed that this 
candidate can accommodate the DAMA, CoGeNT and CRESST experimental results. By 
assuming conservative bounds we have demonstrated that this candidate is not ruled out 
by the CDMS, XENON and PICASSO experiments. Assuming a symmetric Dark Matter 
sector, we also determined the associated thermal relic density and provided relevant bounds 
from indirect search experiments. We found that all the experimental results are compatible 
with a Dark Matter particle with mass around 10 GeV and magnetic moment 1.5 x 10~^® 
e cm, which corresponds to a new physics scale Am ~ 10 TeV. 

Acknowledgments. We thank Robert Foot, Aldo Morselli and Kimmo Kainulainen for 
useful discussions. We are particularly grateful to Marco Regis for sharing parts of the code 
used for the data analysis. 

A Appendix 

In this appendix we estimate the effect of the dipole-dipole interaction term in Eq. (6) for 
all the targets considered in our analysis. To this aim, one can define the useful ratio 

"^^ '^^~ dEn \dEnJ ~ J dasj/dEnvdn^ ' ^ ^ 

which measures the main contribution of the total recoil rate for a given nuclear specie 
in a detector. If r > 1 the interaction is largely spin-dependent, while it is mostly spin- 
independent for r < 1. 

For elastic DM, which features S = 0, the ratio r can be computed from Eqs. (27) 
and (28), yielding 

r{En) = lQsi{Enr(^^^) , (39) 

6 \A Ap rrip J 

where for the illustrative purposes of this appendix we approximated here Fsi(-E_r) ~ 
Fsd(-E_r)- In Fig. 6 we show the ratio r as a function of the nuclear recoil energy considering 
a DM mass of 10 GeV, which gives a good combined fit to the existing experiments, and 
50 GeV, for a comparison. One can see that the dipole-dipole term plays a major role 
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for DAMA (Na, I) and PICASSO (F), while being negligible for all the other experiments 
considered, as stated in the main text. 
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Figure 6: The SD to SI rates ratio r defined in Eq. (38) for the two cases of m^ = 10 GeV 
(left) and m^ = 50 GeV (right). We consider here a Maxwellian halo with local dispersion 
velocity 220 kni/s. Supposing the major part of the signal to come from the lower energy 
threshold (or from the second bin in the case of CoGeNT), r can be determined in the figure 
by the point where the dashed line meets the vertical line with the same color, for each 
nuclear element. 

While the dipole-charge interaction already provides a good fit, the effect of the dipole- 
dipole interaction contributes to a better agreement between the experiments. In fact, the 
increase in the cross section in DAMA favors a lower value for the DM magnetic moment A,^, 
and therefore the DAMA region shifts down towards CoGeNT and CRESST. The overlap 
becomes then almost perfect, for a Maxwellian velocity distribution, especially with velocity 
dispersion 220 km/s. On the other hand also the bound by PICASSO lowers, but even this 
enhanced constraint excludes only a minor part of the overlap zone. 
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